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Abstract 

A method is presented using helicity cross sections for calculating neutrino- nucleon interactions. 
The formalism is applied in the calculation of the pion spectra produced by Va and v T beams. The 
masses of the charged leptons are kept throughout the calculations. Cross sections are presented 
in numerous figures where the contributions of the significant form factors are also shown. The 
article describes the steps of the calculation and gives details so that it can be reproduced and 
adapted to the kinematic conditions of the experiments. 
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I. INTRODUCTION 



Neutrino production of resonances is attracting a lot of attention because differential 
cross sections will be measured in a new generation of experiments which will try to verify 
the functional form of the cross sections i.e. the number and the Q 2 dependence of form 
factors. They will also be used as an input to study properties of neutrinos in oscillation 
experiments. The dominant signal at low neutrino energies will be the A-resonance. Many 
of the completed experiments detected the energy and the angle of the produced muon 
which motivated theoretical authors to integrate over the phase space of the decay products, 
thus presenting cross sections Jw, and the integrated cross section l|, 0, y), Q, 0]. 
The comparisons of the calculated cross sections (differential or integrated) are consistent 
with the data, but we must confess the error bars are large so that there is a significant 
spread on the experimental points. Theoretical calculations of the pion spectra are also 
available 

sag. 

The accuracy will improve in the new experiments and explicit distributions on the 
energy spectrum of the pions produced in the decays will become available. This motivates 
us to calculate the pion spectrum from the diagram in Fig. 1 by keeping the A-propagator 
and without integrating over the whole phase space of the pion. We shall present 
the calculation in detail so that the interested reader can reproduce and use our results. 
We will also make available a code for our calculation which can be used by experimentalists. 



The method that we adopt decomposes the leptonic tensor into helicity components 
and uses helicity cross sections for the scattering of the W ± or Z° bosons on the nucleons. 
This method has been found to be useful ^, S] and was adopted recently in the coherent 
pion production by neutrinos 10J. This way the calculation, whose algebra is long and 
tedious, simplifies. We decided to present results for free protons and neutrons in order 
to show their main features and separate them from nuclear target effects. We also take 
the opportunity to mention and correct a mistake in the pion spectrum that appears in an 
earlier article on which one of us (EAP) was a co-author 3]. 



In addition to describing the formalism we use it for the calculation of the differential 
and integrated cross sections. We calculate the energy spectra of produced pions by 
and v T beams. We keep the masses of the charged leptons throughout the calculations or 
set them equal to zero in order to see the changes brought about in the spectra. We also 
show in many figures the contributions of the important form factors and their interferences 
explicitly. Section [IT] describes the method and includes detailed formulas for the cross 
sections. The functional dependence of the form factors are included in section IHIl where 
they are used in order to calculate the results in figures 2-10. A summary of the results and 
of the improvements that have taken place over the past few years are included in the last 
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section. 

II. THE METHOD 

The process we consider is in general 

v t $) N(p) - r(k')TZ(p A ) - r(k') n(p n ) N(p') (2.1) 
with 1Z being a resonance with spin 3/2. It is convenient to use variables in the rest frame 





FIG. 1: s-channel and u-channel diagrams for the process vip — > £n + p 
of the nucleon 

q = k-k\ Q 2 = -q 2 , W 2 = p 2 A , v = E-E' 
inspired from the kinematics of deep inelastic scattering. The leptonic tensor is 



L =4 



kpk v + k v k^ + g^k.k itfivapk k^ — L^e^ e h 



h,h' 



(2.2) 



(2.3) 



which can be decomposed in terms of the polarizations of the exchanged current. When we 
keep the mass of the muon or tau lepton, there are polarizations for the spin-1 and zero 
states. In the laboratory frame we introduce the basis vectors 



V2 
1 



(0,1,2,0) 



£ £= (0,l,-i,0) 



£ 



(\q\,0,0, go) 



(2.4) 



for helicities and the scalar component 



(2.5) 
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We warn the reader that various definitions occur in the published articles which differ 
from each other. For instance, the above notation is slightly different from than in ref. 10 . 
Numerous articles in the early studies of resonance production and recently [ll|, 0] define 
a set of polarizations in the rest frame of the resonance because it simplifies calculations. 
The above set is convenient with the first three polarizations being present even when the 
leptons are massless and the longitudinal component appearing for massive leptons. It is 
also a complete and orthonormal set of polarizations. 



The coefficients L h i h are obtained by inverting Eq. 12.31 



T — C c-V* cr V 

J-'h'h — >~-\],v^h' t h- 



(2.6) 



When we average over the azimuthal angles of the produced hadrons, only the diagonal 
elements of the density matrix, as well as the £0 interference term survive in the cross 
section. They were calculated in ref. [13] and we give them again for completion: 



J RR 



J LL 



-00 



Q 1 



(2E-u 



2 [Q 2 (2E 



m, 



Tfi 

(2E-u-\q\) 2 - " 



v) 



wm 



212 



Q 2 W 

2ml\Q\2E 



2v{2E-v + \$\) 
2v(2E-v- \q\) + m 

2(Q 2 + ™1) 



v) — wm, 



Q 2 \q\ 



(2.7) 



All matrix elements are positive in the physical region, which becomes evident when the 

is taken into account. For the propagator of the spin- 
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rain 



m. 



kinematic condition Q 
3/2 resonance we introduce the Rarita-Schwinger propagator in free space [13] 



V E-u 



Guv(pa) 



Pa 



Pa 



- Ml + iM A T A 
$a + M a 

- M 2 A + iM A Y A 



9uu 

C 



2 1 



3 Ml 



Pa^Pav 



1 1 



3M/ 



{PAulu ~ PAulu) 



(2.8) 



which is sufficient for the present work. The mass of M A = 1232 MeV and the width 
= 120 MeV will be taken as constants. But in some articles it is a function of the 

( PjrjW) 



120 MeV will be taken as constants 
invariant mass T A = T A0 ( ) with p n (W) 



Several authors studied the modifications of the propagator in nuclear matter, which will 
be useful when we consider nuclear corrections. 
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The matrix element for the entire process includes the coupling of the WpA-vertex given 
in terms of form factors and the npA coupling ig A p%. The matrix element is 



M h = V3 ig A p» u{p' 



G^(p + q) gT + -G u „(p - p w ) d 



u(p) e x (q, h). 



(2.9) 



For our process there are two propagators; one in the s-channel for A ++ , as shown in figure 
1, and one in the w-channel for A . The arguments of the propagator are (p + q) 2 and 
(p ~ Pn) 2 , respectively. The argument of the polarization is q^ and h denotes helicity. In 
this preprint we include, for the calculation and the curves, only the s-channel pole which 
resonates. When we submit the article for publication both terms will be included. 

The coupling at the WpA-vertex are included in d vX which will be discussed below. The 
up A coupling is taken from Appendix Al of ref. [I4J]: g A = 15.3 GeV -1 . The factor \/3 
originates from the isospin relation 



(A ++ lK + b) = v / 3(A + |K 3 |p) 



(2.10) 



where the right hand side is related to the electromagnetic form factor, whose numerical 
value was determined in early experiments and are used in many articles. Electroproduction 
data have been used as an input for neutrino reactions 1 and this convention still survives. 
The factor 1/3 in front of the w-channel pole comes from the Clebsch-Gordan coefficients. 



The WpA-vertex contains vector and axial form factor included in the function 

d vX 
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cj 



M 



N 



+g uX Cf 



qq Ml 



Ml 



7s - q 



riV fW r<v 

3 7 A + + 65 ~ A 



M 



N 



M 



A' 



Ml 



V 



7s 

(2.11^ 



The vector form factors were determined 14j using electroproduction data. Among the axial 
form factor the most important are C^(q 2 ) and C§(q 2 ) and for this reason we omitted the 
other two axial form factors. All form factors will be given explicitly in the next section. 
The square of the matrix element for the s-channel pole is 



M h '*M 



9l 



2 (p 2 A - M 2 A ) 2 + MlY\ 



Tr[ p^A + M A )G' llu d» x e x (h)(p' + M N ) 

e* x/ (h')d v ' x 'G', ul ^ A + M A )p^(f + M N )}{2.12) 



The factor 1/2 comes from averaging over initial spins of the target. The helicity cross 
sections and interference terms for the scattering of the current on a proton target are 



1 To our knowledge the factor v3 was introduced first by Schreiner and von Hippel 16 1 and has become 
traditional to keep it in recent calculations. 
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defined as 

%r<y, Q 2 ) = „ M F M '*M h (2.13) 

We now have all the ingredients for calculating helicity cross sections for the processes 
W + p — > 7£ ++ — > 7r + p. The calculation is straight forward since it involves a two-body phase 
space and a trace. It is long because of the many 7-matrices occurring in the propagator 



and the WpA- vertex. The trace calculation was done using FEYNCALC 15] . 



Finally we can include the lepton variables and present the triple differential cross section 
da G 2 v 



dE w dQ 2 dv 4tt 2 AE 2 



V, 



ad 



\2 



da s da L da R da 1 or da m 



dE^ dE n dE^ dE n dE n 

(2.14) 

This formula includes the muon mass contained in the L^h functions. In the limit m /t = 
it reduces to the known result [8j. We note that the formalism simplifies the calculations 
because the leptonic part was incorporated as an overall factor. We also note that there is 
only one interference term because the other interference terms vanish when we average 
over the azimuthal angle of the produced hadrons. 



III. NUMERICAL ESTIMATES 



Besides the form factors we have now all the ingredients for calculating the pion spectrum. 
The vector form factors have been studied in earlier paper determining their Q 2 -dependence 
from electroproduction data [ijj]. It has been established that they are modified dipoles 

r v (n 2 \ = C s( Q ) 1 to 1 \ 

sWJ (1 + Q 2 /M 2 ) 2 1 + Q 2 /(4M 2 ) [6 - l) 

with C^(Q 2 = 0) = 1.95 and My = 0.84 GeV. The dominance of the magnetic dipole gives 
the relation 

Cl(Q 2 ) = -CUQ 2 )^, C^ = 0. (3.2) 

The other two terms C\ = C\ were set to zero. These form factors were determined 
by electroproduction data where it was shown that they reproduce the measured helicity 



amplitudes [141 ]. 



The axial couplings were obtained from the decay rate of the A-resonance and the repro- 
duction of neutrino data 14 



;i + g 2 /MD 2 i + Q 2 /(3M 

with C 5 A (Q 2 = 0) = 1.2 and M A = 1.05 GeV. PCAC gave us 



Q A (Q 2 ) = ^(Q 2 ) 7 ^ L T- (3-4) 
Q 2 + mi 
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Cq(Q 2 ) is the pseudoscalar form factor with its contribution to the cross section being 
proportional to the square of the lepton mass. 

The triple differential cross section dE d q2 dW shown in figure [2] for neutrino energy of 
1.0 GeV, E n = 300 MeV, Q 2 = 0.2,0.5 and 0.8 GeV 2 . The curves show a A-peak which is 
a sensitive function of Q 2 . This is expected since for large values of Q 2 there is the large 
decrease of the form factors. Integrating over W one obtains the double differential cross 
section of figure [31 Again the cross section decrease with increasing E w because the process 
runs out of phase space. 

Finally more interesting is the dependence on the pion energy when all other vari- 
ables are integrated. Figure H] shows the contribution of the important form factors, 
the term C§ dominates with the next contribution coming from Clf ■ The interference 
between C% and C\ is destructive. The interference between vector and axial form 
factors is constructive for neutrinos and destructive for anti-neutrinos. In figure [5] we set 
the muon mass equal to zero in order to see the effect of neglecting the mass. We re- 
peat the calculations for higher neutrino energies E v = 1, 2, and 5 GeV shown in figures [6J1H1 

The mass of the charged lepton influences the pion spectra shown in Figures H] and El 
Mass effects are much more prominent for v T beams, where the threshold effect is dominant. 
In Fig. [9] we show the pion spectrum for E Vt = 5 GeV. 

A new feature is the change in the significance of the various form factors. The induced 
pseudoscalar form factor C$ is more important relative to C§. The integrated cross section 
is smaller because it reaches its asymptotic value at a much higher energy. In Fig. [10] we 
show the integrated cross section as a function of E Vt . In all the figures, we have included 
the spectra for anti-neutrinos, which are obtained by changing the sign of the vector ® axial 
interference terms. In closing, we remark that the pion spectra show relevant new features, 
and will be important in deciphering the importance and the functional dependence of the 
form factors. 

For comparison with other articles one must keep in mind that we did not include 
nuclear effects from the target. We decided to use protons or neutrons as free targets in 
order to study the significance of the various form factors. We may include nuclear target 
effects later on. 

Several articles calculated and presented the pion spectrum and we comment on them. 
An early article j^, where one of us is a co-author, presented in figures 8-16 pion energy 
spectra with a different shape because the phase space was treated incorrectly. The 
discrepancy was noticed and corrected in figures 3 and 4 of ref. jgj]. The same discrepancy 
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has been pointed out when the spectrum was calculated in [7||. Between the previous 
two articles, mentioned above, and the present article there is a difference in the method 
of calculation. The earlier articles calculated the production of the delta resonance and 
then folded its decay into a pion and a nucleon. This method however requires knowledge 
of the density matrix elements as described in Eqs. (1.3) and (1.4) of ref. [ic|], because 
resonances in various polarization states produce pions with different energies. In the 
present article we calculate the entire process with the A-resonance in the intermediate 
state. The same procedure is advocated in a recent article 0]. Thus the pion energy spec- 
trum on a Hydrogen target is a rather interesting quantity being sensitive to the form factors. 



IV. SUMMARY 



Neutrino interactions are reaching an age of theoretical maturity and will come to 
be compared with the new generation of experiments. For this reason we investigated 
the neutrino cross sections in the energy range of the A-resonance. To this end we have 
rewritten the neutrino-nucleon cross section in terms of helicity of cross sections of the 
on nucleons. For the sake of brevity, we have not included neutral current reactions, which 
will be presented in the future. Our plan is to present these results in a code and explicit 
publications. 



Besides the formalism, we have made the following improvements: 

1. We included the charged lepton mass. The reader who wishes to see effects originating 
form the mass can set them equal to zero in Eqs. 12.71 and include them in the phase 
space of the two-body cross section given in Eq. 12. 131 Such a comparison was presented 
in Figs. H] and [51 

2. We use a running width for the resonance as described after Eq. 12. 131 This was also 



included in ref. 14] 



We study the significance of the pseudo-scalar form factor C$ for muon and tau 
neutrino-induced reactions. We have also presented contributions from the various 
form factors. 



With this article, we hope to clarify several questions presented by colleagues who are 
planning and carrying out the experiments. There are several other quantities that need 
to be calculated, and for just this reason, we have developed a flexible formulation which 
can be adapted to new situations which may arise. Finally, as previously mentioned, we are 
preparing a CODE which will cover new demands that may come up in the future. 
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E v = 1 GeV '- 




FIG. 3: Double differential cross section for a 1 GeV za, interacting with a proton for different E. 
This is calculated from Figl^ after W integration. 
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FIG. 4: Pion energy spectrum for an incoming of energy 0.7 GeV. Contributions from several 
from factors are shown separately. Here, for example, C^&Cq means in this case we have put 
only these two form factors finite. Hence it contains their interference term as well. By Un and 
we denote the contributions of all the form factors to cross sections for these particles. Here muon 
mass effects are taken into account. For a constructive interference is observed, while for it 
is destructive. 
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100 200 300 400 500 600 700 

EJMeV] 

FIG. 5: Same as Fig. ^but with outgoing muon mass neglected. In this limit Cg does not contribute. 
These cross sections are a bit enhanced compared to Fig. [^] due to more phase space in this limit. 
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FIG. 6: Pion energy spectrum for an incoming of energy 1 GeV. Notations are similar to Fig. R[ 
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FIG. 7: Pion energy spectrum for an incoming of energy 2 GeV. Notations are similar to Fig. R[ 



15 



2.50~38 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



2e-38 



©1.5e-38 
o 



LU 
"O 


T3 



1e-38 



5e-39 



'lj C 3 &C 4 'W 



E v =5GeV 




200 400 600 800 
EJMeV] 



1000 1200 1400 



FIG. 8: Pion energy spectrum for an incoming of energy 5 GeV. Notations are similar to Fig. \Q 



16 



1.4e-39 



1.2e-39 



> 
o 



LU 
"O 


"O 



1e-39 



8e-40 



6e-40 



4e-40 



2e-40 



i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i i i i i i i | i i i 



E v = 5 GeV 




200 400 600 800 1000 1200 1400 

EJMeV] 



FIG. 9: Pion energy spectrum for an incoming v T of energy 5 GeV. Notations are similar to Fig. \Q 
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E v JGeV] 

FIG. 10: Variation of v T total cross section with neutrino energy. 
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